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A comparative study of two predictive SO(IO) models , namely the BPW model (proposed by Babu, 
Pati and Wilczek) and the AB model (proposed by Albright and Barr) is done based on their predictions 
regarding CP and flavor violations. There is a significant difference in the structure of the fermion mass- 
matrices in the two models (which are hierarchical for the BPW case and lop-sided for the AB model) 
which gives rise to different CP and flavor violating effects. We include both SM and SUSY contributions 
to these processes. Assuming flavor universality of SUSY breaking parameters at a messenger scale 
M* > Mgut, it has been shown that renormalization group based post-GUT physics gives rise to large 
CP and flavor violations. While these effects were calculated for the BPW model recently, this is the flrst 
time (to our knowledge) that post-GUT contributions have been included for the AB model. The values 
of AmK,ex, Ams^ and S{Bd J/tpKs) are found, in both models, to be close to SM predictions, in 
good agreement with data. Both models predict that S{Bd (l>Ks) should lie in the range -1-0.65-0.74, 
close to the SM prediction and that the EDM of the neutron ~ {few x 10~^®)e-cm, which should be 
observed in upcoming experiments. The lepton sector brings out marked differences between the two 
models. It is found that Br{fj, — > 67) in the AB model is generically much larger than that in the 
BPW model, being consistent with the experimental limit only with a rather heavy SUSY spectrum with 
(mo, mi/2) ~ (1000, 1000) GeV. The BPW model, on the other hand, is consistent with the SUSY 
spectrum being as light as {mo, mi/2) ~ (600, 300) GeV. Another distinction arises in the prediction 
for the EDM of the electron. In the AB model de should lie in the range 10"^^ - 10~^** e-cm, and should 
be observed by forthcoming experiments. The BPW model gives de to be typically 100 times lower than 
that in the AB case. Thus the two models can be distinguished based on their predictions regarding CP 
and flavor violating processes, and can be tested in future experiments. 




Models. 
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1 Introduction 



Grand unified theories [1-3] have found much success in explaining (a) the quantum numbers 
of the members in a family, (b) quantization of electric charge and (c) the meeting of the 
gauge couplings at a scale ~ 2 x 10^^ GeV in the context of supersymmetry [4,5]. In 
particular, it has been argued [6] that the features of (d) neutrino oscillations [7, 8] , (e) 
the likely need for baryogenesis via leptogenesis [9,10], and (f) the success of certain mass 
relations like mi, ~ and m(z/T-)Dirac ~ nT-t at the unification scale, suggest that the effective 
symmetry near the string/GUT scale in 4D should possess the symmetry SU(4)-color [2]. 
Thus, it should be either SO(IO) [11] or minimally G(224) = SU(2)l x SU(2)r x SU(AY [2]. 
(For a detailed review of the advantages and successes of G(224)/SO(10) symmetry, see 
e.g. [6].) 

In recent years, several models based on supcrsymmetric SO(IO) GUT have emerged [12]. 
Two promising candidates have been proposed which have much similarity in their Higgs 
structure and yet important differences in the pattern of fermion mass-matrices. One is 
by Albright and Barr (AB) [13] and the other by Babu, Pati and Wilczek (BPW) [14]. 
Both models use low-dimensional Higgs multiplets (like 45ij, Wh, ^^h and IOh) to break 
SO(IO) and generate fermion masses (see remarks later) as opposed to large-dimensional 
ones (like 126,126,210 and possibly 120). Both of these models work extremely well 
in making predictions regarding the masses of quarks and leptons, the CKM elements 
and neutrino masses and their mixings in good accord with observations. Nevertheless 
there is a significant difference between these two models in the structure of their fermion 
mass matrices. In the BPW-model, the elements of the fermion mass-matrices (con- 
strained by a U(l)-fiavor symmetry [6, 15, 16]) are consistently family-hierarchical with 
"33" > "23" ~ "32" > "22" > "12" ~ "21" > "11" etc. By contrast, in the AB-model, the fermion 
mass-matrices arc lopsided with "23" ~ "33" in the down quark mass- matrix and "32" ~ "33" 
in the charged lepton matrix. (The exact structure of the fermion mass-matrices will be 
presented in Sec. 2.) This difference in the structure of the mass matrices leads to two char- 
acteristically different explanations for the largeness of the — 1/^- oscillation angle in the 
two models. For the BPW model, both charged lepton and neutrino sectors give moderately 
large contributions to this mixing which, as they show, naturally add to give a nearly max- 
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imal sin^ 29^^_^^, while simultaneously giving small Vcb as desired. The largeness of 0^^^^^, 
together with the smallness of Vcb (in the BPW model) turns out in fact to be a consequence 
of (a) the group theory of SO(10)/G(224) in the context of the minimal Higgs system, and 
(b) the hierarchical pattern of the mass-matrices. For the lopsided AB model, on the other 
hand, the large (maximal) — Vj- oscillation angle comes almost entirely from the charged 
lepton sector which has a "32" element comparable to the "33" . 

The original work of Babu, Pati and Wilczek, treated the entries in the mass matrices to 
be real for simplicity, thereby ignoring CP non-conservation. It was successfully extended 
to include CP violation by allowing for phases in the mass matrices by Babu, Pati and the 
author in Ref. [16]. 

The purpose of this paper is to do a comparative study between certain testable predic- 
tions of the AB model versus those of the BPW model allowing for the extension of the latter 
as in Ref. [16]. We find that while both models give similar predictions regarding fermion 
masses and mixings, they can be sharply distinguished by lepton flavor violation, especially 
by the rate of ^ 67 and the edm of the electron. 

We work in a scenario as in Rcfs. [16] and [17], in which flavor-universal soft SUSY 
breaking is transmitted to the sparticles at a messenger-scale M*, with VIgut < M* < 
^string as in a mSUGRA model [18]. Following the general analysis in Ref. [19] it was 
pointed out in Refs. [16] and [17] that in a SUSY-GUT model with a high messenger scale as 
above, post-GUT physics involving RG running from M* — > Mqut leads to dominant flavor 
and CP violating effects. In the literature, however, post-GUT contribution has invariably 
been omitted, except for Refs. [16] and [17], where it has been included only for the BPW 
model. Lepton flavor violation in the AB model has been studied so far by many authors 
by including the contribution arising only through the RH neutrinos [20], without, however, 
the inclusion of post-GUT contributions. I therefore make a comparative study of the BPW 
and the AB models by including the contributions arising from both post-GUT physics, as 
well as those from the RH neutrinos through RG running below the GUT scale. For the sake 
of comparison and completeness, we will include the results obtained in Refs. [16] and [17] 
which deal with CP and flavor violation in the BPW model. 

To calculate the branching ratio of lepton flavor violating processes we include contri- 
butions from three different sources: (i) the sfermion mass-insertions, S^^ arising from 
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renormalization group (RG) running from M* to Mqut ~ 2 x 10^^ GeV, (ii) the mass- 
insertions arising from RG running from M-gut to the right handed neutrino mass 
scales Mij. , and (iii) the chirahty-flipping mass-insertions 5)^^ arising from A— terms that 
are induced solely through RG running from M* to M^j/t involving SO(IO) or G(224) gaug- 
inos in the loop. 

It was found in Ref. [17], that for the BPW- model, contributions to the rate of — > 67 
from sources (i) and (iii) associated with post-GUT physics, were typically much larger than 
that from source (ii) associated with the RH neutrinos. For the AB-model, we find that the 
RH neutrino contribution is strongly enhanced compared to that in the BPW model; as a 
result all the three contributions to the amplitude of /i — 67 are comparable. Including all 
three contributions, we find that for most of the SUSY parameter space, the branching ratio 
for — > 67 calculated in the AB-model is much larger than that in the BPW model and 
is in fact excluded by the experimental upper bound unless (mo, mi/2) ~ 1 TeV. Thus one 
main result of this paper is that, with all three sources of lepton flavor violation included, 
the process /i can provide a clear distinction between the BPW and the AB models. 

We also examine CP violation as well as flavor violation in the quark sector, including that 
reflected by electric dipole moments, in the AB model, and compare it with the corresponding 
results for the BPW- model, obtained in [16]. 

In the following section the patterns of the fermion mass matrices for the BPW and the 
AB models are presented. 

2 A brief description of the BPW and the AB models 

The Babu-Pati-Wilczek (BPW) model 

The Dirac mass matrices of the sectors u, d, I and u proposed in Ref. [14] in the context 
of SO(IO) or G(224)-symmetry have the following structure: 
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(J — 3e 
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r?' - 3e' 

V + 3e' C22 ^-3e 
+ 3e 1 

These matrices are defined in the gauge basis and are multiphed by '^l on left and 
on right. For instance, the row and column indices of are given by (■UL,Ci,fi) 
and (mr, Cr, tR) respectively. These matrices have a hierarchical structure which can be 
attributed to a presumed U(l)-flavor symmetry (see e.g. [6, 16]), so that in magnitudes 1 
^ a Tj ^ C,22 ^ ^22^ > ^' ■ Following the constraints of SO(IO) and the U(l)-flavor 
symmetry, such a pattern of mass-matrices can be obtained using a minimal Higgs system 
consisting of 45h, 16h, 16h, IOh and a singlet S of SO (10)^, which lead to effective couphngs 
of the form [6, 16]: 

>CYuk = /i33l63l63l0,f + [/I23l62l63l0,j(5/M) 

+ a23l62l63l0^^(45^j/M')(5/M)^' + ^23l62l63l6^(16^^/M")(-S/M)«] 

+ [/i2216216210h(^/M)' + (?22l62l62l6^(16j//M") (5/714")"+^] 

+ [^i2l6il62l6^(16^^/M")(^/M)^+2 ^ a,2l6il62l0^^(45ij/M')(5/Mf +2] . (2) 

The powers of (S/M) are determined by flavor-charge assignments (see Refs. [6] and [16]). 
The mass scales M', M" and M are of order Mstring or (possibly) of order Mqut [21]- 
Depending on whether M'{M") ~ Mqut or Mstring (see [21]), the exponent p[q) is either one 
or zero [22]. The VEVs of (45^^) (which is along B - L), {16h) = {IQh) (along {urh)) and 
(S) are of the GUT-scale, while those of (IOh) and (16^) are of the electroweak scale [14,23]. 
The combination 10h-45h effectively acts hke a 120 which is antisymmetric in family space 
and is along B—L. The hierarchical pattern is determined by the suppression of the couplings 
by appropriate powers of Mgut/{M, M'orM"). The entry "1" in the matrices arises from 

^Both the BPW and the AB models beaj: similarities in the choice of the Higgs system, yet there are significant differences 
in the mass matrices. See text for details. 
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the dominant I63I63IOH term. The entries e and e' arising from the 16il6jlOH45H terms, 
are proportional to B — L and are antisymmetric in family space. Thus (e, e')^ — 3(e, e') as 
q ^ I. The parameter a comes from the I62I63IOH term and contributes equally to the 
up and down sectors, whereas fj = rj — a, arising from 16216316^16h operator, contributes 
only to the down and charged lepton sectors. Similarly, (22 arises from the I62I62IOH term 
while C22 g^^s contributions from both I62I62IOH and 16216216^16h operators. Finally, 
1]', which is present only in the down and charged lepton sectors, gets a contribution from 
16i16216^16h terms in the Yukawa Lagrangian (see Eq. (0)). 

The right-handed neutrino masses arise from the effective couplings of the form [24] : 

-CMaj = fijl6d6jlGHl6H/M (3) 

where the /j/s include appropriate powers of {S)/M. The hierarchical form of the Majorana 
mass-matrix for the RH neutrinos is [14]: 



X z 
y 
z y 1 



Mr (4) 



Following flavor charge assignments (see [6]), we have 1 ^ y ^ z ^ x. We expect 
M,t < M < Mpi where M,^ ^ 4 x 10^^ GeV and thus M ^ 10^^ GeV (1/2-2). The 
magnitude of Mr can now be estimated by putting ^ 1, {'^Qh) ~ 2 x 10^^ GeV and 
M ^ (1/2 - 2) lO^s GeV [6,14]. This yields: Mr = h^ljQuf/M ^ (4 x 10^^ GeV) (1/2-2). 

Thus the Majorana masses of the RH neutrinos are given by [6, 14]: 

M3 ^ Mr f« 4 X 10^^ GeV (1/2-2), 

M2 ^ ly^lMa ^ 10^2 GeV(l/2-2), (5) 
Ml ^ [x-z^lMa ~ (l/4-2)10-^M3 
~4 X 10^° GeV(l/8-4). 

Note that both the RH neutrinos as well as the light neutrinos have hierarchical masses. 
In the BPW model of Ref. [14], the parameters a, 77, e etc. were chosen to be real. Setting 
= (2^2 = 0, and with mf^^^ = 174 GeV, m,(m,) = 1.37 GeV, m,(l GeV) = 110 - 116 
MeV, m„(l GeV) = 6 MeV, and the observed masses of e, /x, and r as inputs, for this CP 
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(6) 



conserving case the following fit for the parameters was obtained in Ref. [14] : 

a ^0.110, r/^ 0.151, e ^ -0.095, |?7'| ^ 4.4 x lO'^, 
e' fa 2 X 10-^ M'i fa mt{Mx) ~ 100 GeV, M^^ ~ rnr{Mx) ~ 1.1 GeV. 

These output parameters remain stable to within 10% corresponding to small variations 

(< 10%) in the input parameters of rrit, rric, mg, and m„. These in turn lead to the following 
predictions for the quarks and light neutrinos [14], [6]: 

mbirrib) ^ (4.7-4.9) GeV, 

^ mius) ^ (1/24 eV)(l/2-2), 



Vr 



cb 



rj-e 



rric 


— 1 


mt 


1 a—e \ 



0.044, 



lose 



1/2 



+ 



(7) 



< fa 10.437 + (0.378 ±0.03)1 (for ^ 

Thus, sin2 2^°^';, fa 0.993, 



1/6), 



Vus 



0.20, 



/TO^ ~ 0.07, 

V rric ' 

mail GeV) fa 8 MeV. 

To allow for CP violation, this framework can be extended to include phases for the 
parameters in Ref. [16]. Remarkably enough, it was found that there exists a class of fits 
within the SO(10)/G(224) framework, which correctly describes not only (a) fermion masses, 
(b) CKM mixings and (c) neutrino oscillations [6, 14], but also (d) the observed CP and 
flavor violations in the K and B systems (see Ref. [16] for the predictions in this regard). A 
representative of this class of fits (to be called fit A) is given by [16]: 

(7 = 0.109 - 0.012i, 77 = 0.122 - 0.0464i, e = -0.103, V = 2.4 x lO'^, 

' ' (8) 

e' = 2.35 x 10-V(69°), = 9.8 x lO-^e-^^^^^")^ ^_y^o^ _y^0) ^ (IqO^ 1.1) GeV. 

In this particular fit C22 ^^^^ ^^^^ of economy in parameters. However, 

allowing for (22 ~ (V3)(C22) would still yield the desired results. Because of the success of 
this class of fits in describing correctly all four features (a), (b), (c) and (d) mentioned above 
- which is a non-trivial feature by itself - we will use fit A as a representative to obtain the 
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sfermion mass-insertion parameters S^Il j^j^^, {S^Il)^^'^ ^lr,rl lepton sector and 

thereby the predictions of the BPW model and its enxtension (Ref. [16]) for lepton flavor 
violation. 

The fcrmion mass matrices M^, and M/ are diagonalized at the GUT scale 2 x 10^^ 
GeV by bi-unitary transformations: 

The approximate analytic expressions for the matrices Xf j^^ can be found in [16]. The 
corresponding expressions for X\^ j^ can be obtained by letting (e, e')— > — 3(e, e'). For our 
calculations, the mass-matrices have been diagonalized numerically. 



The Albright-Barr Model 

The Dirac mass matrices of the u, d, 1 and u sectors are given by [13]: 



M„ 










e/3 
-e/3 1 



Mu; Ma 





~5 
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-e/3 



a + e/3 
1 



M 



D 



(10) 
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— e 
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Mi = 
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-e 


Md 
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1 








a + e 1 





These matrices are defined with the convention that the left-handed fermions multiply 
them from the right, and the left handed antifermions from the left. The AB model involves 
a multitude of Higgs multiplets to generate fermion masses and mixings including a 45h, 



two pairs of 16h + 16h, two pairs of IOh and several singlets of SO(IO). The "1" entry 
in the mass matrices arises from the dominant I63I63IOH operator. The e entry arises 
from operators of the form 16216310h45h (as in the BPW model). Since (45h) (x B — L, 
the e entry is antisymmetric, and brings in a factor of 1/3 in the quark sector. The a 
term comes from the operator 16216316h16^ by integrating out the 10s of SO (10). (Note 
that the two 16s of Higgs, 16h and 16^, are distinct). The 16'jj breaks the electroweak 
symmetry but does not participate in the GUT scale breaking of SO (10). The resulting 



operator is 5(162)10(163)(5(16'h))(1(16h)), where the 5,10 and 1 C SU(5). Thus the a 
contributes "lopsidedly" to the 1 and d matrices. The entries S and 6' arise from the operators 
16il6jl6Hl6H, hke the a and contribute only to the 1 and d matrices. Finally, fj, which 

enters the u and u Dirac mass matrices, is of order 10^''"^ and arises from higher dimensional 
operators. The Majorana mass matrix for the right-handed neutrinos in the AB model is 
taken to have the following form: 





2 ~2 


—berj 


at] 




Mr = 


—hefj 




—e 






afj 


— e 


1 





with Ar = 2.5 X 10^^ GeV. The parameters a, b and c are of order one to give the LMA 
solution for neutrino oscillations. Given below is a fit to the parameters a, e, 6 etc. which 
gives the values of the fermion masses and the CKM elements in very good agreement with 
observations [25,26]: 

a =1.78, e = 0.145, 5 = 8.6 x 10-^ 5' = 7.9 x lO'^, 

(12) 

(j) = 126°, f} = 8x 10"^ {Mu, Md) ~ (113, 1) GeV. 
In the next section, we turn to lepton flavor violation. 

3 The Three Sources of Lepton Flavor Violation 

As in Refs. [16] and [17], we assume that flavor-universal soft SUSY-breaking is transmitted 
to the SM-scctor at a messenger scale M*, where Mqut < M* < Mstring- This may naturally 
be realized e.g. in models of mSUGRA [18], or gaugino-mediation [27] or in a class of 
anomalous U(l) D-term SUSY breaking models [28,29]. With the assumption of extreme 
universality as in CMSSM, supersymmetry introduces five parameters at the scale M*: 

mo,mi/2,Ao,taii.(3 and sgn{ii). 

For most purposes, we will adopt this restricted version of SUSY breaking with the added 
restriction that = at M* [27]. However, we will not insist on strict Higgs-squark-slepton 
mass universality. Even though we have fiavor preservation at M*, fiavor violating scalar 
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(mass)^-transitions arise in the model through RG running from M* to the EW scale. As 
described below, we thereby have three sources of lepton flavor violation [16,17]. 

(1) RG Running of Scalar Masses from M* to Mqut- 

With family universality at the scale M*, all sleptons have the mass mo at this scale and 
the scalar (mass)^ matrices are diagonal. Due to flavor dependent Yukawa couplings, with 
ht = hfy = hT-{= h^^) being the largest, RG running from M* to Mqut renders the third 
family lighter than the first two (see e.g. [19]) by the amount: 

Ami = ^^L = ^^f. = = A ^ (t6^^^' ln{M*/MouT) ■ (13) 

The factor 30^12 for the case of G(224). The slepton (mass)^ matrix thus has the form 
= diag( m^, m^, m^ - A). As mentioned earlier, the spin-1/2 lepton mass matrix is 
diagonalized at the GUT scale by the matrices ^. Applying the same transformation 
to the slepton (mass)^ matrix (which is defined in the gauge basis), i.e. by evaluating 
X''1{M.^°^)ll and similarly for L— >R, the transformed slepton (mass)^ matrix is no longer 
diagonal. The presence of these off-diagonal elements (at the GUT-scale) given by: 

i^iL,RRh = {xl^iy^'f^)xQ lm\ (14) 

induces flavor violating transitions — > Here m^- denotes an average slepton mass and 
the hat signifies GUT-scale values. Note that while the (mass) ^-shifts given in Eq. (|13p are 
the same for the BPW and the AB models, the mass insertions 6ll,rr would be different for 
the two models since the matrices Xj^ ^ are different. As mentioned earlier, the approximate 
analytic expressions for the matrices X'l ^ for the BPW-model can be found in [16]. The 
corresponding expressions for X^j^ can be obtained by letting (e, e')^ — 3(e, e'), though we 
use the exact numerical results in our calculations. 

(2) RG Running of the A— parameters from M* to Mgut- 

Even if y4o = at the scale M* (as we assume for concreteness, see also [27]). RG running 
from M* to Mqut induces A— parameters at Mqut, invoving the SO(10)/G(224) gauginos; 
these yield chirality flipping transitions (/}^^ — > V^j^). If we let Mie^ ~ ^ioh ~ ^gut, 
following the general analysis given in [19], the induced A— parameter-matrix for the BPW 
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model is given by (see [17] for details): 



(^Lij)BPW 



M* 



GUT 



-285e' + 90V 

285e' + 90V 90(23-27(22 -285e + 90r/ - 27(t Xjj(15) 
285e + 90r/ - 27a 63 

where Z = (^j^^htgfQMx- The coefficients Y' t) sums of the Casimirs of the 

SO(IO) representations of the chiral superfields involved in the diagrams. For the case of 
G(224), we need to use the substitutions: (f , f , f )^(f , f , f )• The X^ ^ are defined in 
Eq. (|T4|) . The A-term contribution is directly proportional to the SO(IO) gaugino mass Mx 
and thus to mi/2- For approximate analytic expressions of X^^, see Refs. [16] and [17]. 
For the Albright-Barr model, the induced A— matrix for the leptons is given by: 

905 905'e'^ 

906 -95e Xj^ (16) 

905'e^<^ 90a + 95e 63 



lr)ab 



Z ln{ 



M* 



M, 



GUT 



(^l_r)ab is transformed to the SUSY basis by multiplying it with the matrices that diag- 
onalize the lepton mass matrix i.e. X'j^ as in Eq. p5|) . The chirality flipping transition 
angles are defined as : 



"tan (3 m 



(17) 



I ' I 

(3) RG Running of scalar masses from Mqut to the RH neutrino mass scales: 

We work in a basis in which the charged lepton Yukawa matrix Y; and are diagonal 
at the GUT scale. The off-diagonal elements in the Dirac neutrino mass matrix Y^v in this 
basis give rise to lepton flavor violating off-diagonal components in the left handed slepton 
mass matrix through the RG running of the scalar masses from Mqut to the RH neutrino 
mass scales M/j. [30]. The RH neutrinos decouple below M/j.. (For RGEs for MSSM with 
RH neutrinos see e.g. Ref. [31]). In the leading log approximation, the off-diagonal elements 
in the left-handed slepton (mass) ^-matrix, thus arising, are given by: 



\RHN _ (3"^o + ^o) 



J2{Y^MY*),, ln{ 



M, 



GUT, 



The superscript RHN denotes the contribution due to the presence of the RH neutrinos. For 
the case of the AB-model, in the above expression, {YN)ik{Y^)jk {XN)kj{XN)ki because of 
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the definition of tlie mass-matrices. Tlie masses M^. of RH neutrinos are determined from 
Eqs. (0) and for the BPW and AB models respectively. The total LL contribution, 
including post-GUT contribution (Eq. (HH)) and the RH neutrino contribution (Eq. (fTHj) ). 
is thus: 

(sLiir = (^L).. + (^L)r'' (19) 

We will see in the next section that this contribution to yU — > 67 is very different in the two 
models (noted in part in Ref. [33]) and provides a way to distinguish the two models. We find 
that this contribution in the AB model is a factor of ~ 25 — 35 larger in the the amplitude 
than that in the BPW model, and this difference arises entirely due to the structure of the 
mass matrices. We also find that this difference in the mass matrices, also gives rise to large 
differences in the edm of the electron between the two models. 

We now present some results on lepton flavor violation. In the following section we will 
turn to CP violation, and see how the two models compare. 



4 Results on Lepton Flavor Violation 

The decay rates for the lepton flavor violating processes U 1^7 (i > j) are given by: 

m^iti) = ^m\'+\^n (20) 

Here A-'j^ is the amplitude for (/j)^ (^i)^7 decay, while A-'^ = y4mp((/j)^ — ^ The 
amplitudes ^4-}*^ are evaluated in the mass insertion approximation using the (5^^)™, 
and S^Lniii calculated as above. The general expressions for the amplitudes A-'l ^ in one loop 
can be found in e.g. Refs. [31] and [34]. We include the contributions from both chargino 
and neutralino loops with or without the /i— term. 

In Table 1 we give the branching ratio of the process /i ^ 67 and the individual contri- 
butions from the sources 5^*^, ^irjh and {^ij}^^^ (see Eqs. (jlH), (fTTjl and (fTK|) ) evaluated 
in the SO(10)-BPW model, with some sample choices of (mo, mi/2)- For these calculations, 
to be concrete, we set = 1, i.e. M* ^ SMqut, tan/5 = 10, Ao{ at M*) = 

and /i > 0. In the BPW model, for concreteness, the RH neutrino masses are taken to be 
Mr^ = 10^° GeV, Mr^ = 10^2 GeV and Mr^ = 5 x 10^^ GeV (see Eq. 1^). For the masses of 
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the right-handed neutrinos in the AB model, we set M/j^ = 7.5 x 10^ GeV, M/j^ = 7.5 x 10^ 
GeV and Mr.^ = 2.6 x 10^^ GeV corresponding to a = c = 4 and 6 = 6 in Eq. (fTTj) . (The 
results on the rate of /i — » 67, presented in the following table do not change very much for 
other ((9(1)) values of a, b and c). It should be noted that the corresponding values for the 
G(224)-BPW model are smaller than those for the SO (lO)-BPW model approximately by a 
factor of 4 to 6 in the rate, provided ln{M*/MGUT) is the same in both cases (see comments 
below Eqs. / f73|) and / f73)) ). A pictorial representation of these results is depicted in Figs. 1 
and 2. 
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4.5 X 10-13 


-7.2 X 10-11 


3.7 X 10-1^ 


1.3 X 10-^ 


(100, 250) AB 


-8.5 X 10"" 


1.9 X 10-12 


-6.4 X 10-11 


1.3 X 10-12 


8.0 X 10-® 


(500, 250) BPW 


-1.9 X 10-12 


1.0 X 10-12 


-1.6 X 10-12 


8.5 X 10-1^ 


2.2 X 10-11 


(500, 250) AB 


-1.4 X 10-12 


4.4 X 10-12 


-1.4 X 10-12 


2.9 X 10-12 


2.6 X 10-1° 


(800, 250) BPW 


-3.5 X 10-13 


6.1 X 10-13 


-2.9 X 10-13 


4.9 X 10-1^ 


1.3 X 10^2 


(800, 250) AB 


-2.6 X 10^13 


2.5 X 10-12 


-2.6 X 10-13 


1.7 X 10-12 


1.1 X 10-1° 


(1000, 250) BPW 


-1.5 X 10-13 


4.3 X 10-13 


-1.2 X 10-13 


3.5 X 10-1^ 


8.1 X 10^=^ 


(1000, 250) AB 


-1.1 X 10-13 


1.8 X 10-12 


-1.1 X 10-13 


1.2 X 10-12 


5.9 X 10-11 


(600, 300) BPW 


-1.3 X 10-12 


7.2 X 10-13 


-1.1 X 10-12 


5.9 X 10-1^ 


1.1 X 10^1 


(600, 300) AB 


-9.8 X 10-13 


3.0 X 10-12 


-9.7 X 10-13 


2.0 X 10-12 


1.3 X 10-1° 


(100, 500) BPW 


-5.4 X 10-11 


3.5 X 10-1^ 


-2.8 X 10-11 


2.8 X 10-15 


2.6 X 10-« 


(100, 500) AB 


-4.0 X 10-11 


1.5 X 10-13 


-2.5 X 10-11 


9.7 X 10-1^ 


1.6 X 10-^ 


(500, 500) BPW 


-4.3 X 10-12 


3.1 X 10-13 


-3.3 X 10-12 


2.5 X 10-1^ 


1.9 X 10-1° 


(500, 500) AB 


-3.2 X 10-12 


1.3 X 10-12 


-3.0 X 10-12 


8.6 X 10-13 


7.5 X 10-11 


(1000, 500) BPW 


-4.8 X 10-13 


2.6 X 10-13 


-3.9 X 10-13 


2.1 X 10-1^ 


1.4 X 10^2 


(1000, 500) AB 


-3.5 X 10-13 


1.1 X 10-12 


-3.5 X 10-13 


7.3 X 10-13 


1.6 X 10-11 


(200, 1000) BPW 


-1.3 X 10-11 


8.8 X 10-15 


-7.1 X 10-12 


7.2 X 10-1*^ 


1.6 X 10-9 


(200, 1000) AB 


-9.9 X 10-12 


3.7 X 10-1^ 


-6.4 X 10-12 


2.4 X 10-1^ 


1.0 X 10-9 


(1000, 1000) BPW 


-1.1 X 10-12 


7.7 X 10-1^ 


-8.3 X 10-13 


6.3 X 10-15 


1.2 X 10^1 


(1000, 1000) AB 


-7.9 X 10-13 


3.2 X 10-13 


-7.4 X 10-13 


2.2 X 10-13 


4.7 X 10^2 
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Table 1. Comparison between the AB and the BPW models of the various contributions to the amplitude and of 
the branching ratio for ^ ^ 67 for the case of SO(IO). Each of the entries for the amplitudes should be multiplied 
by a common factor slo- Imaginary parts being small are not shown. Only the cases shown in bold typeface are 
in accord with experimental bounds; the other ones are excluded. The first three columns denote contributions to 
the amplitude from post-GUT physics arising from the regime of M* Mgut (see Eqs. II14^ - I|17|| V where for 
concreteness we have chosen \n{M* /Mgut) = 1- The fifth column denotes the contribution from the right-handed 
neutrinos (RHN). Note that the entries corresponding to the RHN-contribution are much larger in the AB-model 
than those in the BPW-model; this is precisely because the AB-model is lopsided while the BPW model is hierarchical 
(see text). Note that for the BPW model, the post-GUT contribution far dominates over the RHN-contribution while 
for the AB model they are comparable. The last column gives the branching ratio of — * 67 including contributions 
from all four columns. The net result is that the AB model is compatible with the empirical limit on /i ^ 67 only for 
rather heavy SUSY spectrum like (mo, mi/2) ~ (1000, 1000) GeV, whereas the BPW is fully compatible with lighter 
SUSY spectrum like (m^, mi/2) - (600, 300) GeV (see text) for the case of SO(IO), and {ruo, mi/2) ~ (400, 250) 
GeV for G(224). These results are depicted graphically in Figs. 1 and 2. 

Before discussing the features of this table, it is worth noting some distinguishing features 
of the BPW and the AB models. As can be inferred from Eqs. ()15|) and ()16|). for a given mo, 
the post-GUT contribution for both the BPW and the AB models increases with increasing 
mi/2 primarily due to the A-term contribution. It turns out that for mi/2 ~ 300 GeV, this 
contribution becomes so large that Br(yU 67) exceeds the experimental limit, unless one 
chooses nio ~ 1000 GeV, so that the rate is suppressed due to large slepton masses. This 
effect applies to both models. 

For the hierarchical BPW model, however, it turns out that the RHN contribution is 
strongly suppressed both relative to that in the lopsided AB-model; and also relative to the 
post-GUT contributions (see discussion below). As a result the dominant contribution for 
the BPW model comes only from post-GUT physics, which decreases with decreasing mi/2 
for a fixed mo. Such a dependence on mi/2 is not so striking, however, for the AB model 
because in this case, owing to the lopsided structure, the RHN contribution (which is not so 
sensitive to mi/2) is rather important and is comparable to the post-GUT contribution. 

Tables 1 and 2 bring out some very interesting distinctions between the two models: 

(1) The experimental limit on /i — > 67 is given by: Br(yU — > 67) < 1.2 x 10^^^ [35]. This 
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means that for the case of the AB model, with dominant contribution coming not only from 
post-GUT physics but also from the RHN contribution, only rather heavy SUSY spectrum, 
{mo, mi/2) > (1000, 1000) GeV, is allowed. The BPW-model, on the other hand, allows 
for relatively low mi/2 (~ 300 GeV), with moderate to heavy rrio, which can be as low as 
about 600 GeV with mi/2 < 300 GcV. As a result, whereas the AB model is consistent with 
II ^ only for rather heavy sleptons 1200 GeV) and heavy squarks 2.8 TeV), the 
BPW model is fully compatible with much lighter slepton masses ~ 600 GeV, with squarks 
being 800 GeV to 1 TeV. These results hold for the case od SO(IO). For the G(224) case the 
BPW model would be consistent with the experimental limit on the rate of ^ 67 for even 
lighter SUSY spectrum including values of (mo, mi/2) ~ (400, 250) GeV, which corresponds 
to rriq ~ 780 GeV and ~ 440 GeV. 

(2) From the point of view of forthcoming experiments we also note that ^ 67 for the 
BPW case, ought to seen with an improvement in the current limit by a factor of 10-50. For 
the AB case, even with a rather heavy SUSY spectrum ((mo, mi/2) ~ (1000, 1000) GeV), 
/X ^ 67 should be seen with an improvement by a factor of only 3-5. Such experiments are 
being planned at the MEG experiment at PSl [32] 

(3) As has been noted earlier in [33] and more recently in [17], the contribution to Al{ij — 
ej) due to RH neutrinos in the BPW model is approximately proportional to 77 — o" ~ 0.041, 
which is naturally small because the entries 77 and a are of C(l/10) in magnitude due 
to the hierarchical structure. In the AB-model on the other hand, this contribution is 
proportional to cr + 2e/3 ^ 1.8. Thus we expect that in amplitude, the RHN contribution 
in the BPW model is smaller by about a factor of 40 than that in the AB model. This has 
two consequences: 

(a) First, there is a dramatic difference between the two models which becomes especially 
prominent if one drops the post-GUT contribution, that amounts to setting M* — Mqut- 
In this case the contribution to (// — > 67) comes entirely from the RHN contribution. In 
this case the branching ratio of {/i — > 67) in the two models differs by a factor of about 
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(40)2 ~ C(103) as depicted in table 2. 



{nio, mi/2)(GeV) 


Br(/i 67)^^ 




(100, 250) 


1.2 X 10"^^ 


9.7 X 10^1^ 


(800, 250) 


2.1 X 10-" 


1.7 X 10-1^ 


(600, 300) 


2.8 X 10-" 


2.5 X 10-1^ 


(500, 500) 


5.3 X 10-^^ 


4.4 X 10-1^ 


(1000, 1000) 


3.4 X 10-13 


2.8 X 10-^^ 



Table 2. Branching ratio for (/i — ^ 67) based only on the RHN contribution (this corresponds to setting M* = 
Mgut) for the AB and BPW models for different choices of (mo, mi/2)- 

It can be seen from table 2 that with only the RHN contribution (which would be the total 
contribution if M* = Mgut)-, the AB model is consistent with the limit on /i ^ 67 for light 
SUSY spectrum, e.g. for (mo, mi/2) = (100, 250) GeV. A similar analysis for the AB model 
was done in Ref. [26] (including the RHN contribution only), and our results agree with 
those of Ref. [26]. One may expect that for the same value of mi/2, increasing nio would 
result in decreasing the branching ratio. For example, from Eq. (jl8p . one may expect the 
rate for /i — > 67 to be proportional to (m^/m^)^ ~ l/^^o- However, the associated loop 
function (see e.g. Ref. [34]) alters the dependence on (mo, mi/2) drastically; it increases 
with increasing mo, for fixed mi/2. The net result of these two effects is that for the same 
mi/2, a low mo ~ 100 GeV and a high nio ~ 1000 GeV, give nearly the same value of the 
branching ratio for fi ^ cy with the inclusion of only the RH neutrino contribution (see Fig. 
3) . This can also be seen in the results of Ref. [26] which analyzes the AB model. The RHN 
contribution in the case of the BPW model is extremely small because of its hierarchical 
structure, as explained above. 

Of course, in the context of supersymmetry breaking as in mSUGRA or gaugino- 
mediation, we expect M* > Mgut, thus post-GUT contributions should be included at 
least in these cases. With the inclusion of post-GUT physics, as mentioned above, the AB 
model is consistent with the experimental limit on /x — 67, only for very heavy SUSY spec- 
trum with (mo, mi/2) > (1000, 1000) GeV, i.e. mj > 1200 GeV and m,- > 2.8 TeV; whereas 
the BPW model is fully compatible with the empirical limit for significantly lower values of 
(mo, mi/2) ~ (600, 300) GeV, i.e. mj ~ 600 GeV and nig ~ 1 TeV (see table 1). 
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(b) Second, it was shown in Ref. [17] that the P-odd asymmetry parameter for the process 
(//+ ^ e+7) defined as A{ii+ e+7) = {\Al\^ - \Ar\^)/{\Al\^ + \Ar\^) (where \Al\ = 
\A^I\6ll) + A^p{6Lfi) + ^l'' I)) is typically negative for the BPW model except for cases with 
very large mi/2 e.g. (m^, mi/2) = (1000, 1000) or (500, 500) GeV. For the AB-casc, due 
to the large RHN contribution, \Al\ > \Af>\ and therefore the P-odd asymmetry parameter 
A would typically be positive. Thus the determination of A in future experiments can help 
distinguish between the BPW and the AB models. 

For the sake of completeness, we give the branching ratios of the processes r — > //7 and 
T — > 67 calculated in the two models in table 3. 



(mo, mi/2)(GeV) 


AB-model 


BPW-model 




Br(r /X7) 


Br(r 67) 


Br(r ^7) 


Br(r — > 67) 


(100, 250) 


2.9 X 10"^ 


3.8 X 10-^^ 


2.6x10^^ 


1.6x10-" 


(800, 250) 


1.0 X 10-s 


4.5 X 10~ii 


1.6x10-9 


6.8x10-12 


(600, 300) 


1.4 X 10"^ 


6.4 X 10-1^ 


2.1x10-9 


8.4x10-12 


(500. 500) 


2.1 X 10-" 


1.0 X 10-^^ 


3.9x10-1" 


1.8x10-12 


(1000, 1000) 


1.5 X 10-^° 


6.5 X 10-^^ 


2.5x10-1^ 


1.1x10-1^ 



Table 3. Branching ratios for (r /i7) and (r 67) evaluated in the two models for the case of SO (10), for some 
sample choices of (mo, mi/2). We have set tan/3 = 10, /j > and In^^ M^iyj ) ~ ^■ 

From table 3 we see that the predictions for the branching ratios for (r — > /I'j) and 
(r — > ej) in either model are well below the current experimental limits. The process 
(r fi'y) can be probed at BABAR and BELLE or at LHC in the forthcoming experiments; 
(r — > 67) seems to be out of the reach of the upcoming experiments. 

In the following section we turn to CP violation in the two models. 

5 Results on Fermion Masses, CKM Elements and CP Violation 

CP violation in the BPW model was studied in detail in Ref. [16]. We will recapitulate 
some of those results and do a comparative study with the AB model. For any choice of the 
parameters in the mass matrices (ry, a, e etc. for the BPW case, and a, e etc. for the AB 
case), one gets the SO(10)-model based values of pw and rjw, which generically can differ 
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widely from the SM-based phenomenological values. We denote the former by (p']y)BPW,AB 
and {ri{^)BPW,AB and the corresponding contributions from the SM-interactions (based on 
p'^r and ?7(y) by SM'. In our calculations we include both the SM' contribution and the SUSY 
contributions involving the sfermion (mass)^-parameters(5^j^ rrlr) "which are in general CP 
violating. These parameters are completely determined in each of the two models for a given 
choice of flavor preserving SUSY-parameters (i.e. ruo, mi/2, fi, and tan/5; we set ^4^ = 
at M*). Using the fits given in Eqs. (jH)) and (fT2|) . we get the following values for the CKM 
elements and fermion masses using mj(mt) = 167 GeV and 171^(171^) = 1.777 GeV as inputs: 
BPW: 

((Ks, Kfe, \Vubl \Vtd\){<mz))BPW- {0-2250, 0.0412, 0.0037, 0.0086) 
{Pw)bpw = 0.150, {r]'w)BPW = 0.374 

(mfe(mfe), m,(m,)) ^ (4.97, 1.32) GeV (21) 

{nisilGeV), m^) ^ (101, 109) MeV 

(mlilGeV), m^ailGeV), m°) ^ (10.1, 3.7, 0.13) MeV 

AB: 

{{Vus, Vcb, \Vub\, \Vu\){<mz))AB- {0.220, 0.041, 0.0032, 0.0081) 
{Pw)ab = 0.148, {r]'^)AB = 0.309 

(mfe(mb), mc(m,)) ^ (4.97, 1.15) GeV (22) 
(msilGeV), m^) ^ (177, 106) MeV 
(mKlGeV), m°(lGe\/), m°) ^ (3.2, 8.5, 0.56) MeV 
The predictions of both models for the CKM elements are in good agreement with the 
measured values, and (p'^y) and f/^) are close to the SM values in each case. It was remarked 
in Ref. [16] that for the BPW model, the masses of the light fermions (u, d and e) can 
be corrected by allowing for O{10~^ — 10~^) "11" entries in the mass matrices which can 
arise naturally through higher dimensional operators. Such small entries will not alter the 
predictions for the CKM mixings. For the AB model, the masses of the bottom and strange 
quarks have been lowered by the gluino loop contributions from 5.12 GeV and 183 MeV 
to 4.97 GeV and 177 MeV respectively. Thus from Eqs. (PT|) and (j221), we see that both 
models are capable of yielding the gross pattern of fermion masses and especially the CKM 
mixings in good accord with observations; at the same time [p'^y) and fj'yy) are close to the 
phenomenological SM values. 
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We now present some results on CP violation. We include both the SM' and the SUSY con- 
tributions in obtaining the total contributions (denoted by "Tot"). The SUSY contribution 
is calculated using the squark mixing elements, which are completely determined 

in both models for any given choice of the SUSY breaking parameters m^, mi/2, Ao, tan (3 
and sgn{fi). As emphasized earlier, in our calculations, the include contributions from 
both post-GUT physics as well as those coming from RG running in MSSM below the GUT 
scale. (For details, see Ref. [16]). We set Ao = for concreteness, as before. Listed below in 
Table 4 are the results on CP and flavor violations in the K° — K° and Bd° — Bd° systems 
for the two models. For these calculations we set lia{M* / Mqut) = 1- 



{ruo, mi/2)(GeV) 


Am^'^-(GeV) 
Tot ^ SM' 


exiSM') 


ex (Tot) 


ArriB, (GeV) 
Tot ^ SM' 


Tot ^ SM' 


(300, 300) BPW 


2.9 X 10-^^ 


2.8 X 10-3 


2.6 X 10-3 


3.5 X 10-13 


0.73 


(300, 300) AB 


2.8 X 10-15 


2.2 X 10-3 


2.2 X 10-3 


3.1 X 10-13 


0.66 


(600, 300) BPW 


2.9 X 10-^^ 


2.8 X 10-3 


2.0 X 10-3 


3.6 X 10-13 


0.73 


(600, 300) AB 


2.8 X 10-15 


2.2 X 10-3 


1.4 X 10-3 


3.1 X 10-13 


0.66 


(1000, 250) BPW 


2.9 X 10-15 


2.8 X 10-3 


1.4 X 10-3 


3.6 X 10-13 


0.74 


(1000, 250) AB 


2.8 X 10-15 


2.2 X 10-3 


-4.0 X 10-3 


3.13 X 10-13 


0.656 


(1000, 500) BPW 


2.9 X 10-15 


2.83 X 10-3 


2.6 X 10-3 


3.6 X 10-13 


0.73 


(1000, 500) AB 


2.8 X 10-15 


2.2 X 10-3 


2.0 X 10-3 


3.1 X 10-13 


0.66 


(1000, 1000) BPW 


2.9 X 10-15 


2.8 X 10-3 


2.9 X 10-3 


3.5 X 10-13 


0.72 


(1000, 1000) AB 


2.8 X 10-15 


2.2 X 10-3 


2.3 X 10-3 


3.1 X 10-13 


0.66 



Table 4. CP violation in the K° — K° and Bd — Bd systems as predicted in the BPW and the AB models for some 
sample choices of {irio, mi/2) and a generic fit of parameters (see Eq.l|SJ for the BPW case and Eq. 11211 for the AB 
case). The superscript s.d. on Attik denotes the short distance contribution. The predictions in either model are in 
good agreement with experimental data for most of the cases displayed above, especially given the uncertainties in 
the matrix elements (see text). It may be noted that values of S{Bd — > J/ipKs) as high as 0.74 in the AB model, 
and as low as 0.65 in the BPW model, can be achieved by varying the fit. 

In obtaining the entries for the i^-system we have used central values of the matrix element 
Bx and the loop functions rji (see Refs. [36, 37] for definitions and values) characterizing 
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short distance QCD effects - i.e. Bk = 0.86 ± 0.13, fx = 159 MeV, rji = 1.38 ± 0.20, 772 = 
0.57±0.01 and 773 = 0.47±0.04. For the S-system we use the central values of the unquenched 
lattice results: Jb^sJ^, = 215(11) (l^g) (15) MeV and fBs\f^. = 245(10) (t^)(+J) MeV 
[38]. Note that the uncertainties in some of these hadronic parameters are in the range of 
15%; thus the predictions of the two SO(IO) models as well as those of the SM would be 
uncertain at present to the same extent. 

Some points of distinctions and similarities between the two models are listed below. 

(1) First note that the data point (mo, mi/2) = (300, 300) GeV displayed above, though 
consistent with CP violation, gives too large a value for Br(/x 67) for both BPW and 
AB models. All other cases shown in table 4 are consistent with the experimental limit on 
/X — > 67 for the BPW model. For the AB model on the other hand, as may be inferred 
from table 1, the choice (mo, mi/2) = (1000, 1000) GeV is the only case that is consistent 
with the limit on // — > 67 (see table 1). It is to be noted that for this case the squark 
masses are extremely high (~ 2.8 TeV), and therefore, in the AB model, once the ji ^ 
constraint is satisfied, the SUSY contributions are strongly suppressed for all four entities: 
Amii", ex, Ams^ and S{Bd J/ipKs). 

(2) For the BPW model on the other hand, there are good regions of parameter space 
allowed by the limit on the rate of ^ 67 (e.g. (mo, mi/2) = (600, 300) GeV), which are 
also in accord with ex- The SUSY contribution to €k for these cases is sizable (~ 20 — 30%) 
and negative, as desired. 

(3) We have exhibited the case (mo, mi/2) = (1000, 250) GeV to illustrate that this case 
does not work for either model as it gives too low a value for ex in the BPW model, and a 
negative value in the AB model. In this case the SUSY contribution, which is negative, is 
sizable because of the associated loop functions which are increasing functions of (m^^/m?). 

(4) The predictions regarding Am^^:, Ame^ and S{B(i — > J/t/jKs) are very similar in both 
the models, i.e they are both close to the SM value. 

(5) As noted above, there are differences between the predictions of the BPW vs. the AB 
models for for a given {rrio, 1711/2)- With uncertainties in B^ and the SUSY spectrum, 
ex cannot, however, be used at present to choose between the two models, but if {rrio, mi/2) 
get determined (e.g. following SUSY searches at the LHC) and Bk is more precisely known 
through improved lattice calculations, ex can indeed distinguish between the BPW and the 
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AB models, as also between SO(IO) and G(224) models (for details on this see Ref. [16]). 
This distinction can be sharpened especially by searches for fi — > ej. 

(6) Bd — > (pKs, '^rriB/. Including the SM' and SUSY contributions to the decay 
4>Ks, we get the following results for the CP violating asymmetry parameter S{Bd (pKs) 
in the two models: 

BPW : S(Bd (j)Ks) ^ +0.65 - 0.74 . 

^ ^ (23) 

AB : S{Bd (pKs) ^ +0.61 - 0.65 . 



The values displayed above for the AB model are calculated for the fit given in Eq. ()12|1 . 
For variant fits in the AB model, values as high as S(5rf (pKs) ~ 0.7 may be obtained. 
The SUSY contribution to the amplitude for the decay B^ — > (pKs in the BPW model 
is only of order 1%, whereas in the AB model it is nearly 5% for light SUSY spectrum 
{{rrio, mi/2) ~ (300,300) GeV) and about 1% for large (m^, mi/2)(~ (1000,500) GeV). 
The main point to note is that in both models S(i?d <pKs) is positive in sign and close 
to the SM prediction. The current experimental values for the asymmetry parameter are 
S{Bd ^ <pKs) = (+0.50 ± 0.25t^oZ) BaBar] (+0.06 ± 0.33 ± 0.09)belle [39]^. While the 
central values of these two measurements are very different, the errors on them are large. It 
will thus be extremely interesting from the viewpoint of the two frameworks presented here 
to see whether the true value of S{Bd (pKs) will turn out to be close to the SM-prediction 
or not. 

Including SUSY contributions to Bg — Bs mixing coming from 5l^j^ ^.i? rl insertions we 
get: 

BPW : /\msXTot ^ SM') ^ 17.3 ps-\^^Y . ^^^^ 
AB: AmsATot ^SM')^ 16.6 ps~\q^^Y . 

Both predictions are compatible with the present lower limit on Am^^ > lAAps"^ [41]. 

(7) Contribution of the A term to e'^: Direct CP violation in Kl vrvr receives a new 
contribution from the chromomagnetic operator Q~ = {g/16n^){sLO'^'^t"'dR — SRa^'^t"'dL)G'^i^, 
which is induced by the gluino penguin diagram. This contribution is proportional to 

At the time of completing this manuscript, the BELLE group reported a new value of S{Bd 4>Ks) = +0.44 ± 0.27±0.05 
at the 2005 Lepton-Photon Symposium [40] . This value is close to the value reported by BaBar, and enhances the prospect of 
it being close to the SM prediction. 
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Im[{S'[f()2i — (^lr)i2]' which is known in both models (see Eqs. (P3|l and (HH)). Follow- 
ing Refs. [42] and [43], one obtains: 

Re{e'/e)g ^ 91 Bg{ ) (— ) Im[{5UU ' {^inYA (25) 

where Bq is the relevant hadronic matrix element. Model-dependent considerations (allowing 
for m\/m^ corrections) indicate that Bq ~ 1 — 4, and that it is positive [42]. Putting in the 
values of 5^^)i2,2i obtained in each model with (mo, mi/2) = (a) (600, 300) GeV, and (b) 
(1000, 1000) GeV, we get: 

BPW : Re{e'/e)~g ^ +(3.7 x 10-^)(5g/4)(10/ tan/5) Case (a) . 

^ +(4.5 X 10-5)(5G/4)(10/tan/3) Case (b) . 

(26) 

AB : Re{e'/e)g ^ -(3.7 x 10-5)(5G/4)(10/tan/5) Case (a) . 

^ +(4.5 X 10-6)(fiG/4)(10/tan/?) Case (6) . 

Whereas both cases (a) and (b) are consistent with the limit on yU — > 67 for the BPW 
model, only case (b) is in accord with ^ 67 for the AB model. The observed value of 
Re{e'/e)obs is given by Re{e'/e)obs = (17 + 2) x 10~^ [44]. At present the theoretical status 
of SM contribution to Re{e'/e) is rather uncertain. For instance, the results of Ref. [45] 
and [46] based on quenched lattice calculations in the lowest order chiral perturbation theory 
suggest negative central values for Re{e'/e). (To be specific Ref. [45] yields Re{e'/e)sM = 
(—4.0 + 2.3) X 10^^, the errors being statistical only.) On the other hand, using methods 
of partial quenching [47] and staggered fermions, positive values of Re{e'/e) in the range of 
about (3 — 13) X 10^^ are obtained in [48]. In addition, a recent non-lattice calculation based 
on next-to-leading order chiral perturbation theory yields Re{e' / e)sM = (19 + 2lg + 6) x 10^^ 
[49]. The systematic errors in these calculations are at present hard to estimate. The point to 
note here is that the BPW model predicts a relatively large and positive SUSY contribution 
to Re{e'/e), especially for case (a), which can eventually be relevant to a full understanding 
of the value of e'^^, whereas this contribution in the AB model is rather small for both cases. 
Better lattice calculations can hopefully reveal whether a large contribution, as in the BPW 
model, is required or not. 

(8) EDM of the neutron and the electron: RG-induced A-terms of the model gen- 
erate chirality-flipping sfermion mixing terms whose magnitudes and phases are 
predictable in the two models (see Eq. (HZj)), for a given choice of the universal SUSY- 
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parameters {m.o, mi/2, and tan/3). These contribute to the EDM's of the quarks and the 
electron by utilizing dominantly the gluino and the neutralino loops respectively. We will ap- 
proximate the latter by using the bino-loop. These contributions are given by (see e.g. [50]): 

(rf.,4)A., = {-I in - 5/(^)/-(<^Sn ^27) 

^ I I 

The EDM of the neutron is given by dn — |(4ci(i — du). The up sector being purely real 
implies d„ = in the AB model. In table 5 we give the values of d„ and de calculated in the 
two models for moderate and heavy SUSY spectrum and tan P — 10. 



{ruo, mi/2)(GeV) 


AB-model 


BPW-model 




dn (e-cm) 


de (e-cm) 


dn (e-cm) 


de (e-cm) 


I (600, 300) 


4.0 X 10-2^ 


1.6 X 10-27 


1.1x10-26 


1.1x10-29 


II (1000, 500) 


1.4 X 10-26 


5.9 X 10-28 


3.9x10-27 


4.1x10-30 


III (1000, 1000) 


5.7 X 10-27 


1.1 X 10-27 


1.7x10-27 


7.7x10-30 


Expt. upper bound 


6.3 X 10-26 


4.3 X 10-27 


6.3 X 10-26 


4.3 X 10-27 



Table 5. EDMs of neutron and electron calculated in the BPW and the AB models for moderate and heavy SUSY 
spectrum and tan /3 = 10 arising only from the induced A-terms. While all cases are consistent with — > 67 for the 

BPW model, only case III is consistent for the AB model. 

Prom the table above, we see that while both models predict that the EDM of the neutron 
should be seen within an improvement by a factor of 5-10 in the current experimental limit, 
their predictions regarding the EDM of the electron are quite different. While the AB model 
predicts that the EDM of the electron should be observed with an improvement by a factor 
of 5-10 in the current experimental limit, the prediction of the BPW model for the EDM of 
the electron is that it is 2 to 3 orders of magnitude smaller than the current upper bound. 
These predictions are in an extremely interesting range; while future experiments on edm of 
the neutron can provide support for or deny both models, those on the edm of the electron 
can clearly distinguish between the two models. 
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6 Conclusions 



In this paper we did a comparative study of two realistic SO (10) models: the hierarchical 
Babu-Pati-Wilczek (BPW) model and the lop-sided Albright-Barr (AB) model. Both mod- 
els have been shown to successfully describe fermion masses, CKM mixings and neutrino 
oscillations. Here we compared the two models with respect to their predictions regarding 
CP and flavor violations in the quark and Icpton sectors. CP violation is assumed to arise 
primarily through phases in fermion mass matrices (sec e.g. Ref. [16]). For all processes we 
include the SM as well as SUSY contributions. For the SUSY contributions, assuming that 
the SUSY messenger scale M* lies above Mqut as in a mSUGRA model, we include contri- 
butions from both post-GUT physics as well as those arising due to RG running in MSSM 
below the GUT scale. While this has been done before for the BPW model in Refs. [16] 
and [17], this is the first time that flavor and CP violations have been studied in the AB 
model including both post-GUT and sub-GUT physics. This inclusion brings out important 
distinctions between the two models. 

Previous works on lepton flavor violation in the AB model [20] have included only the 
RHN contribution associated with sub-GUT physics. It is important to note, however, that in 
both models the sfermion-transition elements (5^^ and the induced A parameters get 

fully determined for a given choice of soft SUSY-brcaking parameters (m^, mi/2, ^oi tan/3 
and sgn{fi)) and thus both contributions are well determined. Including both contributions, 
we flnd the following similarities and distinctions between the two models. 
Similcirities: 

• Both models are capable of yielding values of the Wolfenstein parameters (p'^, r]'^r) which 
are close to the SM values and simultaneously the right gross pattern for fermion masses, 
CKM elements and neutrino oscillations. For this reason, both models give the values of 
ArriK, AiTT-Ba ^iid S{Bfi — * J/ipKs) that are close to the SM predictions and agree quite well 
with the data. The SUSY contribution to these processes is small (< 3%). 

• For the case of €k, it is found that for the BPW model, the SM' value is larger than the 
observed value by about 20% for central values of Bk and rji, but the SUSY contribution 
is sizable and negative, so that the net value can be in good agreement with the observed 
value for most of the SUSY parameter space. For the AB model, for the choice of input 
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parameters as in Eq. (jl2j) . the SM' value for ex is close to the observed value. For most 
of the soft-SUSY parameter space the AB model also yields ex in good agreement with the 
observed value once one allows for uncertainties in the matrix elements (see table 4). 

• Both models predict that Si^B^ (pKs) should be ~ +0.65 — 0.74, close to the SM 
predictions. 

• The predictions regarding Am^^ are similar and compatible with the experimental limit 
in both models. 

• Both models predict the EDM of the neutron to be {few x 10~^^e — cm) which should be 
observed with an improvement in the current limit by a factor of 5-10. 

Thus a confirmation of these predictions on the edm of the neutron and S{Bii —* (pKs), 
would go well with the two models, but cannot distinguish between them. 
Distinctions: 

•The lepton sector brings in impressive distinction between the two models through lepton 
flavor violation and through the EDM of the electron as noted below. 

• The BPW model gives BR(/i 67) in the range of 10^^^ — 10~^^ for slepton masses 
< 500 GeV with the restriction that mi/2 ~ 300 GeV (see remarks below table 1). Thus it 
predicts that fi —>■ e'-f should be seen in upcoming experiments which will have a sensitivity 
of 10^^'^ — 10^^^ [32]. The contribution to /i — 67 in the AB model is generically much 
larger than that of the BPW model. For it to be consistent with the experimental upper 
bound on BR(/i 67), the AB model would require a rather heavy SUSY spectrum, i.e. 
{mo, mi/2) > (1000, 1000) GeV, i.e. mj > 1200 GeV and m^ > 2.8 TeV. With the constraints 
on (mo, mi/2) as noted above, both models predict that fi ^ e'y should be seen with an 
improvement in the current limit which needs to be a factor of 10-50 for the BPW model 
and a factor of 3-5 for the AB model. 

• An interesting distinction between the AB and the BPW models arises in their predic- 
tions for the EDM of the electron. The AB model give de in the range 10~^^ — lO'^'^e cm 
which is only a factor of 3-10 lower than the current limit. Thus the AB model predicts that 
the EDM of the electron should be seen in forthcoming experiments. The BPW model on 
the other hand predicts a value of de in the range 10~^^ — lO^'^^e cm which is about 100-1000 
times lower than the current limit. 

• In the quark sector, another interesting distinction between the two models comes from 
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e'/e. The BPW model predicts that Re{e'/e)susY ~ +5 x 10-4(Sg/4)(10/ tan/3). Thus the 
BPW model predicts that SUSY will give rise to a significant positive contribution to e'/ e, as- 
suming fi^ is positive [42]. The AB model gives Re (e'/e) ~ -5xlO~^(5G'/4)(10/tan/5). 
Thus it predicts that the SUSY contribution is ~ 0(1/10) the experimental value and is 
negative. Since the current theoretical status of the SM contribution to Re(e'/e) is uncertain, 
the relevance of these contributions can be assessed only after the associated matrix elements 
are known reliably. 

In conclusion, the Babu-Pati-Wilczek model and the Albright-Barr model have both been 
extremely successful in describing fermion masses and mixings and neutrino oscillations. In 
this note, including all three important sources of flavor violation (two of which have been 
neglected in the past), we have seen that CP and flavor violation can bring out important 
distinctions between the two models, especially through studies of — > 67 and the edm of 
the electron. It will be extremely interesting to see how these two models fare against the 
upcoming experiments on CP and flavor violation. 
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8 Figures 
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Figure 1: Regions in the {nio, ?7ii/2) plane allowed and disallowed by the current experimental limit on 
Br(/i ej)= 1.2 x 10"" as obtained for the BPW model with ln{M*/MGUT) = 1, tan/J = 10 and ;U > 0. 
The points allowed by the limit on Br(/U — > 67) are marked with a box, while the points disallowed by this 
limit are marked with a star. The results include post-GUT and RHN contributions to the rate of /U — > e-y. 
Note that a large region of parameter space is allowed. 
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Figure 2: Regions in the {nio, 1711/2) plane allowed and disallowed by the current experimental hmit on 
Br(/i ej)= 1.2 x IQ-^^ as obtained for the AB model with ln{M*/MGUT) = 1, tan/3 = 10 and /x > 0. 
The points allowed by the hmit on Br(/U 67) are marked with a box, while the points disallowed by this 
limit are marked with a star. The results include post-GUT and RHN contributions to the rate of — > e-y. 
Note that, only a rather heavy SUSY spectrum with {nio, mi/2) ~ (1000, 1000) GeV is allowed by the limit 
on — > 67. This corresponds to a squark mass of ^ 2.8 TeV and a slepton mass of 1200 GeV. 
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AB Model with only RHN contribution, tan/?=5,yu>0 AB Model with only RHN contribution, tan^= 10,;U>0 
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Figure 3: Curves of constant Br(/i — > ej) in the (mo, mi/2) plane with only the right handed neutrino 
contribution for the case of the AB model. The thickest (blue) line corresponds to the experimental limit of 
1.2 X 10^^^, the medium (green) line to Br(/i — ^ ej) ~ 10^^^, and the thinnest (red) one to Br(/i 67) = 
10"^^. A similar analysis was carried out in Ref. [26]. 
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